NON-NEGATIVE GLOBAL WEAK SOLUTIONS FOR A DEGENERATED 
PARABOLIC SYSTEM APPROXIMATING THE TWO-PHASE STOKES 
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Abstract. We establish the existence of non-negative global weak solutions for a strongly 
couple degenerated parabolic system which was obtained as an approximation of the two-phase 
Stokes problem driven solely by capillary forces. Moreover, the system under consideration 
may be viewed as a two-phase generalization of the classical Thin Film equation. 



1. Introduction and the main result 

In this paper we study the following system of one-dimensional degenerated parabolic equa- 
tions 

dtf = -5. (fdlf + I (2/3 + 3/2^) 93 (/ + g)) , 

dtg = -d. [i^fgdlf + - (2m/ + -if 9 + 6/5') dlif + 5) j 

for (t, x) S (0,00) X X, where X := (0, L) for some L > 0. This system models the motion of 
the interfaces of two immiscible fluid layers of width / and respectively. The layer of width 
/ is located on a impermeable horizontal bottom, identified with the line y = 0, and the layer 
of width g is located on top of the first one. The system (1.1a) has been recently derived in [10] 
as a thin film approximation of the two-phase Stokes problem when capillary is the sole driving 
mechanism. The constants R and /i, which are both assumed to be positive, are determined by 
material properties of the fluids and are given by 

:= ^ and R := ^, 

with fif [resp. fig] denoting the viscosity coefficient of the fluid layer of width / [resp. g]. 
Moreover, 7^ [resp. 7^] is the surface tension coefficient at the interface y = f{t,x) [resp. 
y = f{t,x) + g{t,x)]. The system (1.1a) is supplemented by the initial conditions 

/(0) = /o, 5(0) =50 inX, (1.1b) 

whereby /o and 50 a-re assumed to be known non-negative functions, and by no-flux boundary 
conditions 

d.f = d^g = dlf = dlg = 0, x = 0, L. (1.1c) 
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Let US first observe that if one of the fluid layer has constant zero width, then the system 
(1.1a) becomes, up to a scaling factor, the well-known Thin Film equation 



with m = 3. The theory of existence of weak solutions for the Thin Film equation (1.2) is 
well-established nowadays, cf. [1, 2, 3, 4, 5, 20], to mention just some of the most important 
contributions. We emphasize that it has been rigorously proved in [14, 18] (see also [13]) that 
suitably rescaled solutions of the Stokes and of the Hele-Shaw problem converge towards corre- 
sponding solutions of the equation (1.2) with m = 3 and m = I, respectively. Compared to the 
Thin Film equation (1.2), the system (1.1) is much more complex because it exhibits a strong 
coupling as both equations contain highest order derivatives of all the unknowns. There are also 
two sources of degeneracy because both interfaces may vanish on subsets of the interval I. 

It is worth mentioning that there exists also a two-phase generalization corresponding to the 
thin-film equation (1.2) with m = 1, which has been derived in [9] for flows with capillary and 
gravity effects. The resulting system, which has been investigated in [11, 16, 17] in the presence 
of capillary and in [8, 15] for flows driven only by gravity, appears as the thin layer approximation 
of the the two-phase Muskat problem. Compared to (1.1a), the parabolic system obtained in 
[9] has much more structure: there are two energy functionals available and, furthermore, the 
system can be interpreted as a gradient flow for the L2-Wasserstein distance in the space of 
probability measures with finite second moment. There are not many systems of equations 
which enjoy this nice geometric property. We mention that the parabolic-parabolic Keller-Segel 
system which has a mixed L2-Wasserstein gradient flow structure [6]. As far as we know, the 
two-phase generalization of the thin-film equation (1.2) with m ^ {1, 3} has not been discovered 
yet. 

When studying the problem (1.1), one has to rely only on the energy functional 



a fact which forces us to introduce here a weaker notion of solutions than that in [15, 17]. Indeed, 
it is not difficult to see that the functional £ decreases along smooth solutions of (1.1), as we 
have 



dth = dl 



m > 



(1.2) 




d£{f,g) 



I ((1 + R)d^f + Rd^g)d^{dtf) + R{dM + 9))dAdt9) dx 



dt 



Jx 



[ ((1 + R)dlf + Rdlg)dtf + R{dl{f + g))dtg dx 



Jx 




-iiR^ [ g^\dl{f + g)\^dx- [ f fdlf+^{2f + 'ig)dl{f + g) dx 



Jx Jx 




(1.3) 
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Introducing a suitable regularized version of (1.1a), we construct first, for non- negative initial 
data, globally defined Galerkin approximations which are found to converge towards weak so- 
lutions of the approximating systems. On the other hand, the energy functional E may be 
used to obtain estimates for the solutions of the approximating systems and we obtain sufficient 
information to prove that they converge towards weak solutions of the original problem (1.1). 
Though it is a priori not clear whether the weak solutions of the approximating systems are 
non-negative, we prove that the weak solutions of (1-1) have this property. This differs from the 
framework of thin fluid models with capillary effects and insoluble surfactant [7, 12] where the 
approximating regularized problems may be constructed such that starting from non-negative 
initial data the associated weak solutions are also non- negative. 

Given T G (0, oo] , let Qt '■= (0, T) x X. The main result of this paper is the following theorem, 
establishing the existence of global and non- negative weak solutions for the problem (1.1) that 
start from arbitrary non-negative initial data. 

Theorem 1.1. Let fo,go € H^(I) be two non-negative functions. Then, there exists at least a 
weak global solutions {f,g) of problem (1.1) with the following properties: 

(a) /,gGLoo(0,r;//i(X))n(n,e[o,i/2)C([0,r],C"(X))) forallT>0; 

(b) {f,g){0) = (/o,5o) andf>0,g>Oin (0,oo) x I; 

(c) the mass of the fluids is conserved, that is for every t > 

= ll/olUi and \\g{t)\\Li = Wgoh-,; 

(d) defining for every T > the sets 

Vf := {{t, x)eQT ■■ fit, x) > 0)}, Vg := {(t, x) e Qt ■ g{t, x) > 0)}, 
we have d^f,d^g E L2^ioci'Pf ^'Pg) ^''^^ there exists functions jf,jg,jf,g G L2{Qt) with 
jf = {fd'j + f (2/ + 3fg)d'M + 9)), 

CI 6 tTl J (~] g 

jg = g^l^dl{f + g), jf,g = f'/'gdl{f + g) 
and such that 

Hf := f'/^jf, Hg := f^Rg'^^jg + ^f^'^93f,g + \f^'^93f belong to L^iQr), 



and 



fdt^dxdt+ [ Hfd^^dxdt = [ f{T,x)^{T,x)dx- [ M(0,x)dx, (1.4) 
Jqt Ji Jx 

gdt£,dxdt+ / Hgd^Cdxdt= / g{T,x)C{T,x) dx - / go(,{0,x)dx (1.5) 
Jqt Jx Jx 



for all i e C°°(Qt); 
(e) the energy inequality 



£{f{T),g{T))+ I ^R^\j^\'' + ^\jfj'' + \j^\'dxdt<8{fo,go) (1.6) 



is satisfied for almost all T > 
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We emphasize that due to the lack of regularity of the weak solutions (/, g) found in Theorem 
2.1, which is mainly due to the strong coupling of the system (1.1a), we can identify the function 
Hf only in Vg and Hg only in Vj: 

Hf = (^fd^J + I {2f + 3fg) dUf + g)^ l(o,oo)(/) a.e. in Vg, 

Hg = (^^fgdlf + ^{2f,g' + 3fg + 6fg^)d'M + 9?)\o,oo){9) a.e. inVf. 

Particularly, if the test function in (1.4) satisfies additionally suppSa;^ C Vg, then (1.4) 
is exactly the equation one obtains when multiplying the first equation by (1.1a) by ^ and 
intergrating by parts (similarly for (1.5)). 

The outline of the paper is as follows: in Section 2 we regularize the problem and construct 
global weak solutions for the approximating regularized systems. Furthermore, we prove that 
any accumulation point of the set of approximating weak solutions has to be non-negative. Based 
upon the estimates deduced for this family of weak solutions, we prove in Section 3 that certain 
sequences of approximating weak solutions converge, when letting the regularization parameter 
go to zero, towards non-negative weak solutions of the original problem (1.1). 



2. The regularized approximating problems 

In this section we construct a family of regularized systems approximating in the limit the 
original system (1.1a). This is done in such a manner that the energy functional £ still decreases 
along solutions of the regularized system. The advantage of such a construction is twofold: first 
it enables us to find globally defined Galerkin approximations which are shown to converge 
towards weak solutions of the approximating system, and secondly it provides us with sufficient 
information in order to find accumulation points of this family of weak solutions which solve the 
problem (1.1) in the weak sense defined in Theorem 1.1. 

To proceed we define for every e G (0, 1], the Lipschitz continuous function : M — )• M by 
the relation 

Oe(s) := e + max{0, s} for s G M. 
With this notation, we introduce the following modified version of the original problem (1.1a) 

Otfe = -d.. [alifMfe + f {2al{fe) + ?,al{fe)ae{ge)) dUfe + 5e)) , 

< dtgs = -d, (^^4ife)as{ge)d^Je (2.1) 
+1 {2f,4ige) + 3a2(/,)a,(<7e) + 6ae{fs)4ige)) + 5e)) , 

which is more regular than (1.1a) in the sense that the coefficients of the fourth-order derivatives 
are bounded from below by a positive constant depending only on e. The system (2.1) is 
supplemented by the initial and boundary conditions (1.1b) and (1.1c). 

The main result of this section is the following theorem, ensuring the solvability of the regu- 
larized approximating problem (2.1), (1.1b), and (1.1c) for any e G (0, 1]. 
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Theorem 2.1. Let fo,go G H^(Z) be two non-negative functions and e G (0, 1] he fixed. Then, 
there exists at least a couple of functions {fs,ge) having the following regularity 

• fe,9e G ioo(0, T; H\l)) n L2(0, T; ^^(X)) n (n„e[o,i/2)C([0, T], , 

• dtfs,dtgs£L2{0,T;{HHl))'), 
and satisfying 

j\dtfe{t)\m) dt = 1^ {al{fe)dlf, + I (2a3(/,) + 2,al{fe)aM)) dlife + 9e)^ dxdt, 
j\dt9e{t)m) dt = (I (2/ia3(5,) + 3a2(/,)a,(5.) + Qa,{fe)al{ge)) d^fe + 5.) 

+^a^,{fe)a,{g,)d^Je^ dddxdt 

(2.2) 

for all T > and all ^ G L2{0,T; (Z)), whereby (•!•) is the pairing between H^iT) and 
{H\X)y. Moreover, {fe,gem = (/o,5o), 

jje{T)dx=\\f4L, and j^ge{T) dx = \\go\\L, forallT>0, (2.3) 

dxfe{T) = dxge{T) =0 atx = 0,Lfor almost all T > 0, (2.4) 
and the energy inequality 



n 2 

ae{fe)dlfe + -(2ae(/.) + 3a,(<7.))9^(/. + 



f(/,(r),5e(T))+ / a,(/, 
JQt 

+ fiR'! aU9e)\dUfe+9e)\^dxdt + ^ [ ae{fs)a^M)\d^xife+9 

JQt 4 JQt 



dxdt 



I ^ dxdt 



<£{k,9o) (2.5) 
is satisfied for almost all T > 0. 



2.1. Approximations of the weak solutions of (2.1) by Fourier series expansions. 

In the remaining of this section e G (0, 1] is arbitrary but fixed. In a first step we construct 
Galerkin approximations for the weak solution of the problem determined by (2.1), (1.1b) and 
(1.1c). Letting 

(/)0 := Y^l/L and (f)k := \/2/L cos{k7rx/L), k > 1, 

denote the normalized eigenvectors of the operator —d'^ : H'^iT.) — t- ^2(2^) which satisfy zero 
Neumann boundary conditions, it is well-known that any function which belongs to H^{I) can 
be represented in H^{Z) by its trigonometric series. Let thus /o, go be two non-negative functions 
from H^(I). For each fixed n G N, we consider the partial sums 

n n 

fo-=^hk(t>k, 9o ■=^9(ik(t>k 

k=0 k=0 
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of the series expansions for the initial conditions {fo,go), and we seek for continuously differen- 
tiable functions 

n n 

k=0 k=0 

for t > and x G X, which solve (2.1) when testing with functions from the linear subspace 
{(j)Q, . . . , (pn), and which satisfy initially 

A"(0) = /o", 9e{0)=9l 

By construction the functions {f^,g'^) satisfy the boundary conditions (1.1c) and, if we test 
both equations of (2.1) with (pQ, it follows at once that and are constant in time, that is 

F^{t) = foo, GO(t)=5oo, t>0. (2.6) 

Additionally, testing the system (2.1) successively with (pi, . . . it follows that the 2n— tuple 
(F, G) := (Fg , . . . , F", Gl, . . . , G") is the solution of the initial value problem 

(F,G)' = ^(F,G), (F,G)(0) = (/oi,...,/on,<7oi,...,5on), (2.7) 
whereby the function ^' := (^i, ^2) : I^^" — ^ 1^^" is given by 

^i,i(p,9) =y]pfc / al{Ef{p))dl(pkdx(pjdx 
k=i -^0 

+ i?^^(pfe + qk) ^ [aliEfip)) + ^al{Ej{p))ae{^g{q))^ dl<Pkd^cPj dx, 



and 



'^2,jip,q) =-^Pk / a^,iEf{p))ae{Eg{q))dl(pkd^^jdx 

k=l 

+ RY^Pk + <lk) i^ti4(Eg{q)) + -4(Efip))ae(Eg{q)) 



+ 3ae{Ef{p))al{Eg{q)) ) d'^cpkdxCpj dx 



for J G {1, . . . , n} and (p, q) G M^". We used here the shorthand 

n n 

^f{p) = /oo0o + "^pifpi, ^g{q) = goo(Po + X] ii^i- 
1=1 1=1 

for all p,q £ M". Recalling that is a Lipschitz continuous function, we deduce that ^ is 
locally Lipschitz continuous in M^", and therefore the initial value problem (2.7) possesses a 
unique solution (F, G) defined on a maximal time interval [0,T^). In order to prove that the 
solution is global, that is T" = 00 for all n G N, we show that the energy functional £ decreases 
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along (/^, g^). Indeed, since d'^g^{t) G {<j)Q, . . . , for all t G [0, T"), we may use them 

as test functions for (2.1). Integrating by parts, we then find 

^ {£{f^, g^)) = l^iil + R)d,f^ + Rd^g^)d^{dtf:) + R{d,f^ + d^g^,)d^{dtg^,) dx 
= - + R)dlf: + Rdlg^)dtf^ + R{dl{f^ + g-,))dtg-, dx 



fiR' l^4{g^)\dl{f^ + g^fdx 



3R^ 



as{f^)al{g^)\dl{f^ + g^)\'dx 



R 



ae{f7)dtf7 + ^(2a,(/;) + 2,ae{g:))dt{f^ + g^ 



dx 



[21 



for all t G [0,T"), the last equality following similarly (1.3). Particularly, relation (2.8) ensures 
the boundedness of the function (F, G) on [0, T"). This means that the Galerkin approximations 
{fe^de) exist globally in time for all n G N. 



2.2. Convergence of the Galerkin approximations. 

We next identify an accumulation point of the family ((/", ff" ))n) which is shown subsequently 
to be a weak solution of the regularized system in the sense of Theorem 2.1. To this end, let 
T G (0, oo) be an arbitrary constant. Invoking (2.8), we deduce the uniform boundedness^ of 

d.f^.d^g^, inLoo(0,T;L2(X)), (2.9) 

al'\g-,)dlUe + 97), a\l\f^)aM)dlUe + 9^e) ^ UQt). (2.10) 

al'\f:) [ae{f^)dlf: + |(2a,(/;) + 2,ae{g:))dl{f^ + g^)^ in L^iQr). (2.11) 

while from (2.6) we obtain that 

j^f^{t)dx=\\h\\L, and J^g^{t)dx=\\go\\L, for alH > 0. (2.12) 

First, we observe that Poincare's inequality combined with (2.9) and (2.12) imply that 

f^, g^ are bounded in Loo (0,r;i7i (I)). (2.13) 

On the other hand, by construction we know that > e, and we infer from the relations (2.10)- 
(2.12), by using Poincare's inequality again and the uniform boundedness of {f^)n and (<?")« in 
C(Qt), cf. (2.13), that 

f^, g^ are bounded in L2(0,T;F3(X)). (2.14) 

In the next step, we derive uniform bounds for the time derivatives of the Galerkin approx- 
imations. In order to do so, we observe that the relation of (2.1) may be written in a more 



All the bounds in this section are uniform in n £ N. 



J. ESCHER AND B.-V. MATIOC 



concise form as dtf^ = —dxH'^/^, whereby we set 



Is - '^x^^f 

H}'- ■.=al{f^)dlf: + I {2al{f:) + ^al{f:)aM)) dlUe + 9e) 



The relations (2.11) and (2.13) yield now that the sequence {H^j:'"')n is bounded in L2{Qt)- 
Given ^ G H^(I), we define for each n G N the truncation 

n 
fe=0 

Integration by parts then implies that 

m^m)L,\ =\{dtf:{t)\c)L,\ = \{H}'''{t)\d.Qn)L,\ < wH^'^'imacwm 

<\\H^fmLM\\m- 

Consequently, for every t G [0,T], the function dtf^{t) belongs to the dual {H^{Z)) of H^{Z) 
and, integration with respect to time, yields 

\\^tfe\\L2{0,T;(m{X))') ^ II "f^/ " II L2 (Qt) • 

We claim that a similar estimate is valid also for dtg^ . Indeed, the second relation of (2.1) may 
be recast as the equation dtg^ = —dxHg'^, whereby 

H^' :=^«?(/r)«e(5r)5'/r + f i^f^^ig^) + 3alif^)aM) + Qasif^)4{g^)) dlife + 9e) 

+ hl/\f^)aM) + f (2«.(/r) + 3a,(g,")) dUf^ + g^) 

(2.15) 

Gathering (2.10), (2.11), and (2.13), we see that (i?g'") „ is a bounded sequence in L2{Qt)j and 
we finally conclude that 

{dtf^)n, {dtg^)n are bounded in L2(0, T; {H\Z))'). (2.16) 

Using an argument based on the Aubin-Lions Lemma, cf. Corollary 4 in [19], together with 
the continuity of the embeddings 

comv - for a G [0,1/2), 

we conclude from (2.14) and (2.16) that the sequences {f^)n and {g^)n are both relatively 
compact in C([0, T], C"(X)) n L2(0, T; C2+"(X)) for ah a G [0,1/2). Hence, using a diagonal 
procedure, we find functions and g^ and subsequences of (/") and ((7") (not relabeled) such 
that 

f:^h and g^.^g, in C7([0, T], C"(X)) n L2(0, T; C72+"(X)) (2.17) 
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for all a G [0,1/2). Furthermore, let us observe that the relations (2.13), (2.14), and (2.17) 
ensure that the limit functions belong also to /e,<^e G Loo{0,T; {Z)) n ^2(0, T; i7'^(X)) and 
that 

dlf^^dlfe and dlg'^^dlge in L2(Qr) for A; = 1, 2, 3. (2.18) 
Finally, from (2.16) we obtain that dtfe, dtge G ^2(0, T; (iJi(X))'), and 

dtf^ - dtfe, dtg^ ^ dtge in L2{0,T; {H\X))'). (2.19) 

The fact that {fe,ge) can be defined globally follows by using a standard Cantor diagonal 
argument (choosing a sequence T„ 00). 

2.3. Construction of the weak solutions for the regularized system. 

In this last part of Section 2 we prove that the functions {fs,ge) constructed in (2.17) are 
weak solutions of (2.1), (1.1b), and (1.1c), and enjoy all the properties stated in Theorem 
2.1. First, let us observe that the functions fe and ge possess the regularity and integrability 
properties required in Theorem 2.1. Moreover, because /o G H^{I) and due to (2.17), we have 
that fe{0) = fo and ge{0) = go- Furthermore, combining (2.12) and (2.17), it follows that the 
identities (2.3) are satisfied. 

Concerning (2.4), we note that (2.17) guarantees that /"(T) — )■ fe{T) in C^iX) for almost all 
T > 0. Because dxf^{T) = at x = 0, L, the desired claim (2.4) for fe (and similarly for ge) is 
immediate. 

We next prove that the energy estimate (2.5) is satisfied by the functions {fe, ge)- To this end, 
we infer from the relations (2.10), (2.11), (2.17), and (2.18), by using also the Lipschitz continuity 
of the map and after extracting further subsequences of (/") and {g"^) (not relabeled) that 
we have the following weak convergences in L2{Qt) ■ 

4/\g^)d^Me + 97) - al/\ge)dlUe + 9e), 
al/Hf^)ae{g^)dlif^ + 5?) ^ al/Hfe)ae{ge)dlife + ge), 

a\'\f:) {ae{f:)dlf: + ^{2ae{f:) + 3a,(<7r))9^(/," + g^)^ 

^al/\fe) {ae{fe)dlfe+'^{2ae{fe) + ^ae{ge))dl{fe+ge)^ ■ 

Recalhng that /"(t) — ;> fe{i), dei^) — ^ 9e{t) in C^(X) for almost ah t > 0, and that the initial 
data fo,go belong to H^{I), we conclude after integrating (2.8) with respect to time and passing 
to liminf„_!.oo that the desired energy inequality (2.5) is satisfied. 

To finish the proof of Theorem 2.1, we are only left to prove relations (2.2). Let therefore 
^ G L2{0,T; H^{I)) be given, and define for each ?i G N the truncation 

n 

C{t,-) ■.= ^m-)\(t>k)L,4>k, te{o,T). 

Using integrating by parts, we find in a similar way as before that 

j^dtf^{t)C{t)dx = - I^C{t)d.Hy''{t)dx = I^H}'''{t)d,e{t)dx, 
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whence, we have 
fT 



]mce by Le 
together with (2.19), that 



{dtm\e{t))L, dt 



(2.20) 



/o Jo JQt 

for all n > 0. Since by Lebesgue's dominated convergence — )• ^ in L2{0,T; H^[I)) we find. 



^{dtfmcmdt^ r{dtMt)m)dt. (2.21) 

Jo 
Furthermore, (2.10), (2.13), (2.17), and (2.18) ensure that, after extracting further subsequences, 
in L2{Qt), whereby we set 



we have H^'^ 



R 



{fMfe + ^ {2aUfe) + 3aUfe)ae{9e)) dl{h + 9e 



(2.22) 



Letting now n — )■ oo in (2.20), we obtain from (2.21) and (2.22) the first relation of (2.2). On 
the other hand, combining the relations (2.10), (2.11), (2.13), (2.17), and (2.18) we may assume 
H^g in L2{Qt), with given by 



that Hg'"' 



R 



-al{h)ae{ge)dlfe + ^ (2/ia3(5,) + 3a2(/,)a,(g,) + Qae{fe)al{ge)) dHh + Qe). 

Repeating the arguments presented above we conclude that the second identity of (2.2) is also 
satisfied, and the proof of Theorem 2.1 is complete. 

Let us remark that we do not know whether the weak solutions {fe,ge) found in Theorem 
2.1 are non-negative. The next lemma though, together with the convergence results that we 
will provide in the next section yields the non- negativity of the weak solutions of problem (1.1), 
which are found as being accumulation points of the family {{fe, 9e))ee(0,i]y cf. Lemma 2.2 
and Corollary 2.3 below. To this end, we introduce the following notation. Pick a function 
(fi G C°°(M), which is non-negative, has support contained in [—1,0], and satisfies 



(p{x) dx = 1. 

JR 

Moreover, let the function xi : M — t- IR be defined by the relation 

rX />00 

Xi{x) ■= ~ / tpir) drds for x G M, 

Jo Js 

and {x5)s>o be the associated mollifier, that is Xsi^) •= ^Xii^/^) ^oi x € M and 6 > 0. The 
following properties of {x5)5>o play a-n important role in the proof of Lemma 2.2: 

11X5 - max{-id,0}||i^(K) < 6, (2.23) 

< 1. Wxsh^im ^^'^yha.m^ ^^^'i llxsh^m < ^~^\\f'\\L^{R) (2-24) 



IXsWLo 



for aU (5 > 0. 

Lemma 2.2. The functions {f^,g^) found in Theorem 2.1 satisfy 



X^,{fe{T))dx 



for all e £ (0, 1] and all T > 0. 



< CVfe and 



Xy^{9e{T))dx 



< CVTe 



(2.25) 
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Before proving Lemma 2.2 let us draw the conclusion that all accumulation points of the 
family ((/e, 5e))e6(o,i] in C{Qj',M?), with T > 0, are non-negative functions. 

Corollary 2.3. Assume that there exists a sequence {£k)k C (0, 1] with \ and a pair 
{f,g) G C(Qj.,M2) such that 

{f,^,g,^)^{f,g) tnC(QT,M.^). (2.26) 
Then, f and g are both non-negative functions in Qt- 
Proof. In virtue of (2.24), we have 

IIXv?^(/£j -max{-/, 0}||i^(Q^) 

< Wx^ifsk) - X^iDhooiQT) + IIXv^(/) -max{-/,0}||i^(Q^) 

< Wfck - f\\L^(QT) + IIXv^(/) -max{-/,0}||i^(Q^). 

Whence, our assumption (2.26) guarantees the convergence X,/ej:{fek) ~^ nia-xj— /, 0} in C{Qx)- 
Letting now /c — )• oo in the first inequality of (2.25) yields 



j m&x{-f{t),0}dx = 



for all t G [0, T]. This is the desired assertion for /. The proof of the non-negativity of g follows 
similarly. □ 

Proof of Lemma 2.2. Let 5 > be given. Since x'sifei'^)) ^ ^^i^)^ compute that 

I / Xsif^{t))dx = [ X'sif:{mfsit)dx= j dtf^{t)j2{Me{tmk)LAkdx 

n 

H'/''{t)Y,d. iix'sifemMLAk) dx, 

k=0 

relation which is satisfied for all t > 0. The assertions (2.25) are obviously true when T = 0, 
so let us assume that T > 0. Integration the previous identities with respect to time on [0, T] 
shows that 

Xs{f^{T))dx= [ xs{f7mdx+ [ Hy''f2ix5{feMk)L,d,<t^kdxdt. (2.27) 
Jx Jqt 

In order to let n — >• oo in (2.27), we first observe 

n 

J2(^'s(fe)\ML,d,cPk ^ X5ife)d.fe m L2{Qt). (2.28) 

k=0 

Indeed, we have 

n / " \ 

Xs{fe)d.fe - Y.(MMk)L,d,<l^k = Xs{fe)d.fe " Y.^x'5ifeMk)L,d^cl>k 

k=0 \ k=0 / 

n 

k=0 
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and the convergence of the first term to zero foUows by using Lebesgue's dominated convergence 
theorem together with the fact that x'sifei'^)) ^ ^^i^) t > 0. On the other hand, the 

reminding sum is the truncation of the Fourier series of Xs{f^)9xfe ~ Xsife)9xfe a-^d, using 
(2.24), may be estimated as follows 



k=0 



<\\X5{fe)(^xf£ - X5{fs)dxfE\\L2{QT) 



L2{Qt) 



<2||X5(/") - Xi{fe)\\Loc{QT)\\^^f^\\L2{QT) 

+ 2iix^'(/r)iiL(Q,)ii9./r-a./,iii^( 

<25 ^ll¥''lli^(R)ll/" - fe\\\^{QT)\\^xfe\\L2{QT) 



2{Qt) 



+ 2'^ '^\W\\\^{W.)\\dxf£ - dxfe\\\2{QT)^ 

the desired estimate (2.28) being now a consequence of (2.13) and (2.17). 

Thus, letting n — )• oo in (2.27) and taking into account that /^(O) = /o > 0, we obtain the 
following identity for the weak solution of (2.1) found in Theorem 2.1 



X5{fe{T))dx 



Hfx5{fe)dxfedxdt. 



Since = on M \ (—5, 0), Holder's inequality leads us to 



Xs{fe{T))dx] < 



< 



'S<h<0] 

asift 



H'f\x'Kfe)\dxfe\dxdt 



S<fe<0] 



ae{fe)dlfe + |(2a,(/,) + ?,ae{ge))dl{fe + Qe 



dxdt 



4{fe)x's{fe?\dxfe?dxdt. 

■S<fe<0] 

We choose now S := y/e. Recalling that = e on (— oo,0], the energy inequality (2.5) together 
with (2.24) imply that 



X^Me{T))dx 



1/2 



< c 



e'x':^{fer\dM 



-V^<fe<0] 



dxdt 



/ \dxfefdxdt) <CVfe, 
Jqt J 



which is the desired estimate (2.25) for /g. Concerning the second estimate of (2.25), similar 
arguments to those presented above yield that 



X5{ge{T))dx= / mxsi9e)dx9edxdt 
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for all T > and 6 > 0. Writing Hg as the sum of three terms, cf. (2.15), we obtain from 
Holder's inequality and the estimate (2.5) the following inequalities 



X&{9e{T))dx 



< fJ^R 



3R 



+ 



-S<ge<0] 



a'jH9e)x'R9eM9e\ al'\g,)dl{h + 9e) 



dxdt 



■s<ge<o] 

al/\f^)aM)xK9e)\d.9e 



dxdt 



S<9e<0] 



aV\f:) ( ae{f:)dtf: + ^ (2a,(/;) + ^aM)) d^f^ + 



R 



|3/' 



dxdt 



< C 



1/2 



al{g,){a,{fe) + ae{ge))x'5{9e?\d^9e? 
and, when 5 = y^, we arrive at the following estimate 



/ \d^ge\'^ dxdt ] <cVfe. 
Jqt / 



X^e{9e{T))dx 
This proves the lemma. 

3. Existence of weak solutions for the original problem 



□ 



This last section is devoted to the proof of our main result Theorem 1.1. Therefore, we 
collect first some estimates for the family of weak solutions ((/e, S'e))£6{o,i] of the approximating 
problems (2.1), (1.1b), and (1.1c). Considering now e G (0,1] as a parameter, we deduce from 
(2.2), (2.3), and (2.5) the uniform boundedness of 

d,fe, d,ge in Loo(0,r;L2(X)), (3.1) 

dtfe, dtge mL2(Q,T-{H\X))'), (3.2) 

aT{ge)dl{fe + ge): al/\f,)ae{ge)d'Afe + 9e) m L^iQr), (3.3) 



al/\fe) ( ae{fe)dtfe + ^(2a,(/,) + ?,ae{ge))dt{fe + ge) ) m L2(Qt 



R 



(3.4) 

for all T > 0. Recalling also (2.3), the arguments used in the previous section ensure the 
existence of a sequence (efc)fc C (0, 1] with et \ and functions 

/,5 G L2(o,r;/7i(x))n (n„e(o,i/2]C([o,r],c"(x))) , 

having the property that 

fe.^f, 9e,^9 in C([0,r],C°(X)) for aha G [0,1/2); (3.5) 
fe,-f\ 9e,-9 mL2{0,T;H\l)). (3.6) 

(3.7) 
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Particularly, (3.5) implies that for almost every t G [0, T] we have 

d^fe,{t) ^ d.fit), d^ge,{t) ^ d^g{t) inL2(X), (3.8) 

and it follows now directly from (2.5) that f,g£ Loo{0,T; {!)). Mor eover, the convergence 
(3.5) together with the Corollary 2.3 yield the non-negativity of the limits / and g. The latter 
property combined with the relation (2.3) ensure the desired mass conservation property claimed 
by Theorem 1.1 (c). Let us also observe that the claim (b) of Theorem 1.1 is a simple consequence 
of the convergence (3.5) and of the relations /^(O) = /o and (7e(0) = for all e S (0, 1]. 

We next establish the identities (d) of Theorem 1.1. For this let ^ G C°^{Qx) be given and, 
for every e G (0, 1], let {{f^,g2))n be the sequence found in Section 2 to converge towards the 
weak solution {fe,ge) of problem (2.1). Integrating by parts, we then find that 

r [ dtf:(dxdt= [ f^{T,x)C{T,x)dx- [ f^m{0,x)dx- [ r f^dt^dxdt. (3.9) 
Jo Jx Jx Jx JxJo 

Hence, letting n — )■ oo in (3.9), we deduce in virtue of (2.17), (2.19), and of the first identity in 
(2.2) that 

/ Hfd^Cdxdt= I fe{T,x)C(T,x)dx- [ fo{x)aO,x)dx- [ fedt^dxdt. (3.10) 
Jqt Jx Jx Jqt 

Particularly, recalling (3.5), it is easy to see that along the subsequence ((/ej,,5'ej,))fc, the right 
hand side of the equality (3.10) converges towards the corresponding quantities appearing in the 
equation (1.4). 

In order to study the behaviour of the left hand side of the relation (3.10), we define for every 
m G N, m > 1, the open subsets 

:= {(t, x) G (0, T) X X : f{t, x) > 1/m}, 
K ■= ^) e (0, T) X X : g{t, x) > 1/m}, 

of Qt, and observe that Vf = VJ'^^^Pf and Vg = U^^^P™. Let now m > 1 be fixed. Due 
to (3.5), we may find a positive integer Uq with the property that f^^{t,x) > (2m)~^ and 
gek{t,x) > (2m)-^ for ah {t,x) G Vf nV^ and ah k > k^. Thanks to (3.3) and (3.4), the 
sequences (i9^/ei,)fc and (9^5£j,)fc are both bounded in L2(^™ 1^^™)) Siwd, up to the extraction 
of a diagonal subsequence, we may assume that 

dlfe,-dlf, dig,, ^ dig in L^iVf nV^) (3.11) 

for all m > 1. Because {H^)}^ and {a,^^"^ {fe,)!!^)}^ are bounded in L2{Qt)^ cf. (2.3), (3.1), 
and (3.4), we can also presuppose that there exist functions Hf,jf G L2{Qt) such that 

Hf^Hf, {ae,{fe,))-'''^Hf ^jf mL2{QT). (3.12) 

Using the convergences (3.11) and (3.5), we may identify the weak limits in (3.12) in the set 
where / and g are both positive 

Hf = fdlf + f (2/3 + 3fg) dlif + g) ■ ^ . 

Jf = f'/mdlf + §{2f + 3fg)dl{f + g)) " ''^ 
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Moreover, because of \asi^{fef,) — f \ < £fc + Ifet ~ f\ for sll k > 0, we may identify Hf in the 
large set Vg. Indeed, by the dominated convergence theorem ae{^(/ej.) — t- /^^^ in L2{Qt)^ which 
shows, together with (3.12), that Hf = f^^'^jf in L2{Qt)- Summarizing, we have shown that 

Hf=(^fdlf + ^{2f + 3fg)dl{f + g)^l(^o,oo){f) hi L^iVg), (3.13) 

and the desired assertion (1.4) follows now at once. The identity (1.5) is obtained by using 
similar arguments. Indeed, in this case it is possible to identify first the weak limit Hg of (a 
subsequence of) {Hg'')^ in L2{Qt)- On the other hand, because of (3.3), there exist functions 
iff'i/.s ^ L2{Qt) such that 

4i^{9e,)dlUe,+9e,)-jg, ali\f,Ja,,{gs,)dl{fe,+ge,)^jf,g in L2{Qt). (3.14) 

Again, due to (3.5) and (3.11), we identify jg = g^/'^dl[f + g) and jf^g = f^/'^gd^{f + g) in 
L2{Vf r]Vg). Writing Hg'' in a similar manner as in (2.15), the dominated convergence theorem 
shows then 

Hg = ^Rg'/'jg + ^-^f'^jf^g + If^'^jf in L2{QtI 

and therefore 

Hg=(^-fgdlf + ^[2i,g^ + ?,fg + Qfg^)dl{f + g)^ l(o,oo)(5) in L2(P/). 

The assertion (1.5) is now immediate. 

Finally, we collect (3.8), (3.12), and (3.14), and pass to liminf^,_>oo in the energy inequality 
(2.5) to obtain the desired claim (e) of Theorem 1.1. 

Because T was chosen arbitrary, we may again pick a sequence T„ oo and, extracting a 
diagonal sequence of {{fek^9ek))k we may assume that / and g are globally defined and the 
claims of Theorem 1.1 are true for all T > 0. 

References 

[1] F. Bernis. Viscous flows, fourth order nonlinear degenerate parabolic equations and singular elliptic problems. 

In Free boundary problems: theory and applications (Toledo, 1993), volume 323 of Pitman Res. Notes Math. 

Ser., pages 40-56. Longman Sci. Tech., Harlow, 1995. 
[2] F. Bernis and A. Friedman. Higher order nonlinear degenerate parabolic equations. J. Differential Equations, 

83(l):179-206, 1990. 

[3] A. L. Bertozzi and M. Pugh. The lubrication approximation for thin viscous films; the moving contact line 

with a "porous media" cut-off of van der Waals interactions. Nonlinearity, 7(6): 1535-1564, 1994. 
[4] A. L. Bertozzi and M. Pugh. The lubrication approximation for thin viscous films; regularity and long-time 

behavior of weak solutions. Comm. Pure Appl. Math., 49(2):85-123, 1996. 
[5] M. Bertsch, R. Dal Passo, H. Garcke, and G. Griin. The thin viscous flow equation in higher space dimensions. 

Adv. Differential Equations, 3(3):417-440, 1998. 
[6] A. Blanchet and Ph. Laurengot. The parabolic-parabolic Keller-Segel system with critical diffusion as a 

gradient flow in R'', d > 3. preprint arXiv: 1203.3573. 
[7] J. Escher, M. Hillairet, Ph. Laurengot, and C. Walker. Weak solutions to a thin film model with capillary 

effects and insoluble surfactant. Nonlinearity, 25:2423-2441, 2012. 
[8] J. Escher, Ph. Laurengot, and B.-V. Matioc. Existence and stability of weak solutions for a degenerate 

parabolic system modelling two-phase flows in porous media. Ann. Inst. H. Poincare Anal. Non Lineaire, 

28(4):583-598, 2011. 



16 



J. ESCHER AND B.-V. MATIOC 



[9] J. Escher, A.-V. Matioc, and B.-V. Matioc. Modelling and analysis of the Muskat problem for thin fluid 

layers. J. Math. Fluid Mech., 14:267-277, 2012. 
[10] J. Escher, A.-V. Matioc, and B.-V. Matioc. Thin-film approximations of the two-phase Stokes problem. 

Nonlinear Anal., 2012. DOI:10.1016/j.na.2012.07.034. 
[11] J. Escher and B.-V. Matioc. Existence and stability of solutions for a strongly coupled system modelling 

thin fluid films. NoDEA Nonlinear Differential Equations AppL, 2012. DOT 10.1007/s00030-012-0166-l. 
[12[ H. Garcke and S. Wieland. Surfactant spreading on thin viscous films: nonnegative solutions of a coupled 

degenerate system. SIAM J. Math. Anal, 37(6):2025-2048 (electronic), 2006. 
[13[ L. Giacomelli and F. Otto. Rigorous lubrication approximation. Interfaces Free Bound., 5(4):483-529, 2003. 
[14] M. Giinther and G. Prokert. A justification for the thin film approximation of Stokes fiow with surface 

tension. J. Differential Equations, 245(10):2802-2845, 2008. 
[15[ Ph. Laurengot and B.-V. Matioc. A gradient flow approach to a thin film approximation of the Muskat 

problem. Calc. Var. Partial Differential Equations, 2012. DOT 10.1007/s00526-012-0520-5. 
[16[ Ph. Laurengot and B.-V. Matioc. A thin film approximation of the Muskat problem with gravity and capillary 

forces. 2012. preprint. 

[17[ B.-V. Matioc. Non-negative global weak solutions for a degenerate parabolic system modeling thin films 
driven by capillarity. Proc. Roy. Soc. Edinburgh Sect. A, 142 (5):1071-1085, 2012. 

[18[ B.-V. Matioc and G. Prokert. Hele-Shaw fiow in thin threads: A rigorous limit result. Interfaces Free Bound., 
14(2):205-230, 2012. 

[19[ J. Simon. Compact sets in the space L*'(0,T;_B). Ann. Mat. Pura Appl. (4), 146:65-96, 1987. 
[20[ A. Tudorascu. Lubrication approxiiuation for thin viscous films: asymptotic behavior of nonnegative solu- 
tions. Comm. Partial Differential Equations, 32(7-9):1147-1172, 2007. 

Institut fur Angewandte Mathematik, Leibniz Universitat Hannover, Welfengarten 1, 30167 
Hannover, Germany. 

E-mail address: escher@ifam.uni-hannover.de 

Institut fur Mathematik, Universitat Wien, Nordbergstrasse 15, 1090 Wien, Osterreich 
E-mail address: bogdan-vasile .matiocOunivie . ac . at 



